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LHCb has collected the world’s largest sample of charmed hadrons. 
This sample is used to search for direct and indirect CP violation in 
charm. Recent and updated measurements from several decay modes are 
presented. 
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1 Introduction 


Studies of CP violation are a primary goal of the LHCb Experiment at CERN. CP 
violation in charm decays is expected to be very small in the Standard Model [^[^. 
However, asymmetries at a few times 10“^ within the Standard Model cannot be 
excluded |4]-^. A signihcant excess of CP violation with respect to the theoretical 
predictions would be a signature of physics beyond the Standard Model. The study of 
CP violation in singly Cabibbo-suppressed (SCS) charm decays is uniquely sensitive 
to physics beyond the Standard Model, in particular through new contributions to 
AC* = 1 strong penguin and chromomagnetic dipole operators [^. The analysis of 
the following decay channels allows us several opportunities to probe both integrated 
and localized CP violation in different regions of decay phase space. Large CP 
violation has been observed in regions of meson decay phase space not associated to 
a resonance, further motivating localized searches [^. To date, no CP violation has 
been observed in the charm quark system. 

2 Searches for direct CP violation in charm 

2.1 DC 

(s) A 

The measured time-integrated CP asymmetry for a meson {Q) decay to a hnal state 
/ can be given by: 


Araw(Q /) 


N{Q ^ /) - N{Q ^ /) 
N{Q^ f) + N{Q^ fY 


( 1 ) 


where N indicates the number of reconstructed events of a given decay after back¬ 
ground subtraction. We study the decays t , where h = n or K. CP 

violation is only expected to be measurable in the SCS subset of these decays. The 
analysis is performed using promptly produced D mesons from the full 3 fb“^ LHCb 
Run 1 data sample, which consists of 1 fb“^ of proton-proton collision data taken 
in 2011 at y/s = 7 TeV, and 2 fb“^ of proton-proton collision data taken in 2012 at 
y/s = 8 TeV. 

These measured (or “raw”) asymmetries are a sum of several independent asym¬ 
metries: 


= Acp{K^shC + + ^cp/inu ( 2 ) 

where Acp is the physical CP asymmetry, Ap is the production asymmetry of the 
decay particle, Ap, is the detection asymmetry of the charged daughter hadron, and 
Acpiint represents the detection asymmetry between a and a K^, due to the 
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presence of mixing and CP violation in the kaon system and the different interaction 
rates of and in the detector material [^. 

The quantities of interest can then be determined by forming differences that 
result in the cancellation of production and detection asymmetries: 

Acp{D+ ^ ^ KlA) - ^ 07r+) (3) 

and 

AcpiD^ ^ KlK+) = {A^AD^ ^ KlKA - A^Dj ^ ( 4 ) 

- (klrawP+^i^V)-^rawW^0Vr+)), 

where we assume the Standard Model expectation of no CP violation in the Cabibbo- 
favored (CF) D'^ —>■ (j)7r~^ decay mode. 

As a cross check, a double difference of asymmetries largely insensitive to produc¬ 
tion and instrumental asymmetries may also be formed: 

A^^ = HAUDt ^ KlA) - (A^Dt ^ KIkA) (5) 

- {{A^AD+ ^ KlA) - {AUD+ ^ K^KA) 

^ Acp{Dt ^ A^^7r+) + Acp{D+ ^ K^KA- 

In Equations |3]j^ above, the neutral kaon asymmetry is not shown, but is accounted 
for in the hnal results. As the production and detection asymmetries depend on 
the kinematic distributions, weights are assigned to the candidates such that the 
kinematic distributions are equalized. 

The analysis is performed separately for 2011 and 2012 data, and for both LHCb 
magnet polarities; consistent results are obtained. The combined results 

Acp{D+ K^KA = (0.03 ± 0.17(stat) ± 0.14(syst))%, (6) 

Acp{D^ —)■ = (0.38 ± 0.46(stat) ± 0.17(syst))%, (7) 

and 

Acp{D+ + ^cp{Dt K^A) = (0.41 ± 0.49(stat) ± 0.26(syst))%, (8) 

show no indication of CP violation. These are the most precise measurements of 
these quantities. 

2.2 -A Ah- (A^cp) 

A search for a time-integrated CP asymmetry in Ah~ decays is performed 

using the full LHCb Run 1 data sample (3 fb~^). The flavor of the initial state 
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is tagged by the charge of the muon in semileptonic B —)■ decays, where 

5 is a B~ or B^. The measured asymmetry for tagged neutral D mesons to a CP 
conjugate hnal state / is given by: 


^raw(/) — 


N{B D^fi-X') - N{B D^^i+X') 


(9) 


Af(B -» A'') + N(B D»ii+X') 

where N indicates the number of reconstructed events of a given decay after back¬ 
ground subtraction and X' represents the set of undetected hnal state particles from 
the semileptonic B decay. 

The measured asymmetry is a sum of the physical CP asymmetry {Acp{f)), the 
production asymmetry {Ap{B)) and detection asymmetry (Ad(//)): 

^raw(/) = Ac'p(/) + Ap{B) + Ap,{fl). (10) 

An experimentally robust variable may be constructed, AAcp, by taking the differ¬ 
ence of the measured asymmetries in —)■ K^K~ and —)■ tt+tt" decays: 


AAcP = Araw(A'A') - Araw(7r7r), 


( 11 ) 


as the production and the muon detection asymmetries cancel to hrst order. 

Alternatively, for extracting Acp{KK), the detection and production asymme¬ 
tries can be measured using CF B —>■ Zl°(—)■ K~t[^)^~X decays where no CP vio¬ 
lation is expected. An additional detection asymmetry, AD(A'7r), arises due to the 
different interaction rates of the charged kaon with matter. To remove this asym¬ 
metry, the CF control channels and —)■ (for which we 

assume negligible CP violation) are used. In the D~^ — )■ decays, Acp/mt 


(see Section 2.1) is estimated from simulation and subtracted from the measured 
asymmetry. Once AAqp and Acp{KK) are measured, the individual asymmetry 
Acpinn) = AAcp — Acp{KK) can be computed. 

The analysis is performed separately for the 2011 and the 2012 data, and for the 
two magnet polarities; consistent results are obtained. The combined results for the 
asymmetries [^, 

( 12 ) 

(13) 

(14) 


AAcp = (+0.14±0.16(stat) ±0.08(syst))%, 
Acp{KK) = (—0.06 ± 0.15(stat) ± 0.10(syst))%, 
and, with correlation p = 0.28, 

Acp{'K'k) = (—0.20 ± 0.19(stat) ± 0.10(syst))%, 


are compatible with CP conservation. 

The Acp{hh) asymmetries are the most precise measurement of these individual 
asymmetries to date. The precision of AAcp is comparable to the preliminary result 
for AAcp measured using prompt decays reconstructed in 1 fb“^ of LHCb Run 1 


data from 2011 10 
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2.3 ^ TT'^TT TT^ (Energy Test) 

The energy test 0 is an unbinned model-independent statistical method than can 
be used to search for time-integrated CP violation in —)• 7r“7r+7r° decays. This 

method relies on the comparison of and flavor samples and is sensitive to CP 
violation localized in the phase-space of the multi-body hnal state. 

The flavor of the prompt is tagged by the charge of the slow pion in the 
decay D* —)■ For the reconstruction of the two categories of reconstructed 

neutral pions are used: pions for which both hnal state photons are reconstructed 
separately (resolved pions), as well as pions that have higher momentum (typically 
Pt > 2 GeV/c) and thus a smaller opening angle of the two photons (merged pions). 

The energy test is used to assign a p-value for a non-zero CP violation hypothesis. 
In this method, a test statistic T is used to compare the average distances based on 
the metric function 'ijj. The test statistic is dehned as 


j ^ i’ij Wij _ 

rn(n — 1 r) — 


A 




i,j>i 


n[n 


1 ) 


i,j>i 


n(n 


1 ) 


nn’ 


(15) 




and the metric function 'ifjij = is chosen as a Gaussian function 

with a tunable parameter a. T compares the average distances of pairs of events 
belonging to two samples of opposite havor. The normalization factor removes the 
impact of global asymmetries. The distance between two points in phase space is 
given by dij = (m^g “ ''^ 121'"^23 ~ '"^ 231'"^13 ~ '"^ 13)1 where the subscripts {1,2,3} 
indicate the hnal-state particles. 

If there is GP violation, T is expected to be larger than zero. This technique 
calculates a p-value under the hypothesis of GP symmetry by comparing the nominal 
T value observed in data to a distribution of T values obtained from permutation 
samples, where the flavor of the is randomly reassigned to simulate samples with¬ 
out CP violation. The p-value for the no CP violation hypothesis is obtained as the 
fraction of permutation T values greater than the nominal T value in the data. 

We study the data sample of 2 fb“^ collected by LHGb during 2012, and hnd a 
p-value of (2.6 ± 0.5) x 10“^ 
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This result is based on 1000 permutations. The 
data sample has been split according to various criteria to test the stability of the 
results. Analyses of sub-samples with opposite magnet polarity, with different trigger 
configurations, and with hducial selection requirements removing areas of high local 
asymmetry of the tagging soft pion from the D*^ decay all provide consistent results. 
Various checks have been performed to ensure there are no asymmetries arising from 
background events or detector related asymmetries. Varying the metric parameter a 
results in similar p-values at the 10“^ level. 

This analysis provides the world’s best sensitivity from a single experiment to 
local CP violation in P° —)■ decays. 
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2.4 > K^K TT+TT (Triple product asymmetries) 

Multi-body D decays are sensitive to CP violation due to the interference of sev¬ 
eral resonances across the multi-body phase space. In decays, 

triple products of hnal state particle momenta in the rest frame, dehned as 
Ct = Pk+ ■ {Pn+ ^ Ptt-) Ct = Pk- ■ {Pn- ^ Pit+) mesous, respec¬ 

tively, are quantities which are P-odd (and T-odd). The decay rate asymmetries: 


_ r(C'T > 0) - r(C'T < 0) 
^ “ r(C'T > 0) + r(C'T < 0) ’ 


^ _ t{-Ct > 0) - r(-gT < 0) 

^ “ T{-Ct > 0) + T{-Ct < 0) 


(16) 


are thus sensitive to CP violation. However, hnal state interactions signihcantly alter 
the measured asymmetries, and so the difference = ^{At — At) is used to access 

the CP asymmetry of the meson by canceling hnal state interaction ehects. 

The observable is, by dehnition, insensitive to production and detection 

asymmetries, and is thus very robust against systematic nncertainties. Using the fnll 
LHCb Rnn 1 data sample (3 fb“^), the phase space integrated measnrements of the 
asymmetries are 13 : 


At = (-7.18 ± 0.41(stat) ± 0.13(syst))%, 
At = (-7.55 ± 0.41(stat) ± 0.12(syst))%, 


(17) 

(18) 


and 


^ T-odd 
Icp 


= ( 0.18 ± 0.29(stat) ± 0.04(syst))%. 


(19) 


These results show no evidence of global CP violation in these decays, bnt achieve a 
signihcant improvement in precision over previous measurements 
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The same measnrements are also performed in 32 bins of the hve-dimensional 
phase space. The asymmetries are extracted in each bin and a^p’^'^ calculated. A 
test for consistency across the phase space is performed, yielding a p-value of 74%. 
Thus there is no evidence for direct local CP violation in this decay. 

Similarly, binning in the decay time of the P° candidates and performing the 
same test gives sensitivity to indirect CP violation. This yields a p-value of 72%; no 
evidence for indirect CP violation was fonnd. 


3 Searches for indirect CP violation in charm 

3.1 -A h+h- {Ar) 

A measurement of the indirect CP violation in mixing can be performed in 
the stndy of two-body hadronic charm decays to CP eigenstates (P° —)■ K^K~ or 
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Figure 1: Raw CP asymmetry as a function of decay time for —)■ K^K decays 
(left) and of ^ 7i^7i~ decays (right), and corresponding pull plots 15 . 


—)■ 7r’''7r“). It can be evaluated by the asymmetry in the effective lifetimes (r) of 
flavor-tagged decays and expressed by the following equation with the assumption of 
negligible direct CP-violating contribution: 


, t{D^ ^ f) - t{D^ ^ f) 1 ^ . 

Ar = —=-~ - Am y cos 6 — x sm < 

r(D0 ^/) + r(H0 ^/) 2 ^ ^ 


( 20 ) 


where Am is dehned by \q/p\^‘^ ~ IP Am- A measurement of Ap differing signihcantly 
from zero is a manifestation of indirect CP violation as it requires a non-zero value 
for Am or 0 = arg(g/p). 

A measurement of Ay at LHCb is performed using 3 fb“^ of data from the full 
LHCb Run 1 data sample 15 . In this analysis, Zl° decays are obtained from B —)■ 


D^fi Ufj^X decays and the neutral D meson is tagged by the charge of the muon. Ai 
is determined through a fif to the time-dependent raw CP asymmetry A^(t): 


A^Sp{t) - Ay, - Ay-. 


( 21 ) 


The raw CP asymmetry is determined from simultaneous hts to rrih+h- in bins of 
decay time for both Zl° and Zl° samples, as seen in Figure 

Systematic uncertainties include uncertainty in the Zl° mass £t model, the decay 
time resolution, mixing in the beauty quark system, detection and production asym¬ 
metries, and the time-dependent efficiency. The main contributions to the systematic 
uncertainty result from the mistag probability (the probability that the neutral D 
meson was tagged with the wrong flavor) and the mistag asymmetry (differences in 
positive and negative muon mistag probability), which are studied using a control 
sample of —)■ K~tt^ decays. The mistag asymmetry is the dominant contribution 
to the systematic uncertainty. 
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We find: 


AriKK) = (-0.134 ± 0.077 


AriTnr) = (-0.092 ± 0.145 


-0.0347 /O 

)%. 


( 22 ) 

(23) 


where the uncertainties are statistical and systematic, respectively. This measurement 
is consistent with no indirect CP violation in the charm system, and has comparable 
precision to LHCb’s previous result using neutral D mesons from prompt D* decays 
obtained from 1 fb~^ of LHCb Run 1 data from 2011 (Tol. 


3.2 “Wrong-sign” -A Ktt {CP violation via mixing) 

At LHCb, charm mixing parameters are determined by the time-dependent ratio of 
—)■ K^tt~ (“wrong” sign, WS) to —)■ iC“7r+ (“right sign”, RS) decay rates. 

The RS decay rate is dominated by its CF amplitude. The WS rate arises from the 
interfering amplitudes of the doubly Cabibbo-suppressed decay (DCS) and the CF 
decay following —D^ oscillation. 

Assuming no CP violation and small charm mixing parameters (|a:| and \y\ 1), 

this ratio is: 


m 


rd + ‘ ‘ 


r 


r 


(24) 


where x' = xcos6 + ysm6, y' = ycos6 — xsintf, is the ratio of suppressed-to- 
favored decay rates, and 6 is the strong phase difference between the DCS decays and 
the CF decays {AiD^^ K+tt-) /AiD^^ K-7r+) = - 

We can search for CP violation in these decays by comparing the time-dependent 
ratios evaluated separately for and D^. A difference in the Rd parameter for 
and would be a sign of direct CP violation, while a difference in and y' 
parameters would imply indirect CP violation (|g/p|7^1or07^O). The data are fit 
under three hypotheses: CP is conserved, only an indirect CP contribution is allowed, 
and both direct and indirect CP contributions are allowed. 

The full LHCb Run 1 data sample (3 fb”^) is used to perform these measurements. 
The sample is obtained from D*^ —)■ decays, which allow the determination of 

the flavor of the neutral D meson. Assuming CP conservation, the mixing parameters 
are measured to be Rd = (3.568 ± 0.058 ± 0.033) ■ lO-^, y' = (4.8 ± 0.8 ± 0.5) ■ lO-^ 
and = (5.5 ± 4.2 ± 2.6) ■ 10“® where the first uncertainty is statistical and the 
second systematic 17 . In the other scenarios, the direct CP asymmetry is Ad = 
(—0.7 ± 1.9)%, and the magnitude of q/p is determined to be 0.75 < \q/p\ < 1.24 
and 0.67 < \q/p\ < 1.52 at 68.3% and 95.5% confidence level (C.L.), respectively. 
The results of the measurements of mixing parameters, and of the conhdence regions, 
are shown in Figure for the three scenarios. The results are compatible with CP 
conservation and provide the most stringent bounds on the parameter \q/p\ from a 
single experiment. 
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Figure 2: Two-dimensional confidence regions in the y') plane obtained (a) allow¬ 
ing both direct and indirect CP violation, (b) assuming no direct CP violation and 
(c) assuming CP conservation. The solid (dashed) curves in (a) and (b) indicate the 
contours of the mixing parameters associated with {D^) decays. The solid, dashed 
and dotted curves in (c) indicate the contours of CP-averaged mixing parameters at 
68.3%, 95.5% and 99.7% C.L. The best-fit value is shown with a point 
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